Pulse propagation in decorated granular chains: An analytical approach 
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We study pulse propagation in one-dimensional chains of spherical granules decorated with small 
grains placed between large granules. The effect of the small granules can be captured by replacing 
the decorated chains by undecorated chains of large granules of appropriately renormalized mass and 
effective interaction between the large granules. This allows us to obtain simple analytic expressions 
for the pulse propagation properties using a generalization of the binary collision approximation 
introduced in our earlier work [Phys. Rev. E in print (2009); Phys. Rev. E 69, 037601 (2004)]. 



PACS numbers: 46.40.Cd,43.25.+y,45.70.-n,05.65.+b 



I. INTRODUCTION 

It is well known that an initial impulse applied at one 
end of a granular chain in the absence of gaps or precom- 
pression can result in solitary waves propagating through 
the medium P, Apart from addressing fundamen- 
tal problems of pulse dynamics in the presence of highly 
nonlinear interactions, this scenario also has direct prac- 
tical application, e.g.,in designing shock absorbers (EGI j 
sound scramblers [Ell and actuating devices Q- Thus, 
in the recent past pulse propagation in one-dimensional 
(ID) chains of granules has been studied extensively both 
theoretically and experimentally H, [ p. [ToL [i ll ItH . fl3l . [Til . 

HE EE IE EE HE HE HE EE EE a MlT 

Several configurations and paramctrizations of the ID 
chain have been used to study pulse propagation in the 
face of dissipation [IE HE HE HE] and of polydispersity 
in the structure and mass of the granules [E [lj, H3, 
EH, EE EE El EE EE E3- In particular, polydisper- 
sity is frequently introduced in a regular fashion such 
as in tapered chains (TCs), in which the size and/or 
mass of successive granule s sy stematically decreases or 
increases [E EE El M El EE EE E3 ■ Polydispersityis 
also introduced by distributing masses randomly p8l . E^] , 
by "decorating" chains with small massesplaced reg- 
ularly or randomly among larger masses [l7l . |2?| . and 
by optimizing grain distribution [30l | for particular pur- 
poses. In a recent study [H E3 it was shown that the 
shock absorbing properties of a granular chain can be en- 
hanced quite dramatically by decorating the chain with 
smaller granules. To our knowledge, however, most stud- 
ies on decorated chains have been numerical and hence 



restricted to specific parameter values. One exception is 
the early work of Nesterenko and Laizardi [3l| based on 
a continuum approximation, which is difficult to imple- 
ment for chains with rapid spatial variations such as are 
inherently present in decorated chains. It is clearly desir- 
able to devise new theoretical means to obtain analytic 
results for these systems. This is the goal of the present 
work. 

In a recent paper [32| we analyzed pulse propagation in 
a variety of undecorated TCs. Using a binary collision ap- 
proximation, we presented analytical expressions for var- 
ious quantities that characterize propagation of the pulse 
or pulse front in ID TCs. However, an essential clement 
of the success of this approximation is that a given pair 
of granules collides only once, that is that there be no re- 
bounds. This approach is therefore not directly applica- 
ble to decorated chains, where one (or more) small/light 
granules between two large/heavy ones will rattle back 
and forth. A similar conclusion about the inapplicabil- 
ity of a two-body collision approximation (different from 
ours) to decorated chains was also reached in Ref. [27j . 
On the other hand, the simplicity offered by the binary 
collision approximation is too tempting to discard it al- 
together. Instead, we have succeeded in extending this 
theoretical scheme to decorated chains. 

Thus, in this work we present such a scheme that al- 
lows us to make use of the binary collision approximation 
for decorated chains in which single small/light granules 
are placed between large/heavy ones. The case of more 
than one small granule placed between large ones is not 
yet part of our scheme, nor is it in most decorated con- 
figurations that have been considered (numerically) in 
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the literature. Such a generalization is however entirely 
feasible. Our approach is based on a mapping of the dec- 
orated chain onto an effective undecorated chain. The 
effect of small decorating granules is included through a 
modified interaction between the bigger granules. This 
also requires the renormalization of the mass of the big- 
ger granules. We present a systematic way to carry out 
this renormalization. The binary approximation is then 
used on the effective chain, which allows us to obtain 
analytical results for the quantities that characterize the 
pulse front propagation in these chains. We compare the 
analytical results with exact numerical results for two dif- 
ferent decorated chains: simple decorated chains, where 
the sizes of the bigger particles are all the same, and dec- 
orated linear forward TCs, where the sizes of the bigger 
particles decrease linearly along the chain. 

In Sec. |TT] we introduce the general granular chain 
model in terms of rescalcd (dimensionless) variables. In 
Sec. IIIII we introduce decorated chains and present an 
effective description in terms of the renormalized inter- 
action and masses. We then apply the binary collision 
approximation in Sec. IIVI and obtain analytical results 
for pulse properties. The analytic results are compared 
to numerical integration results in Sec. [V] In Sec. IVII 
we provide a final summary and some comments about 
ongoing work. 



II. THE MODEL 

We consider chains of granules all made of the same 
material of density p. When neighboring granules col- 
lide, they repel each other according to the power law 
potential 



V = -r' k \y k - 2/fc+iT 



(1) 



Here y k is the displacement of granule k from its position 
at the beginning of the collision, and a is a constant deter- 
mined by Young's modulus and Poisson's ratio [H, 34 1 . 
The exponent n is 5/2 for spheres (Hertz potential) 34 1. 
We have defined 



R'k 



fe+1 



1/2 



(2) 



where R' k is the principal radius of curvature of the sur- 
face of granule k at the point of contact. When the gran- 
ules do not overlap, the potential is zero. The equation 
of motion for the fcth granule is 



d 2 Vk 
dr 2 



ar' k -i(Vk-i - Vk) n %fc-i - Vk) 
ar , k (y k -y k+1 ) n - 1 9(y k -y k+1 ), (3) 



where M k = (4/3)7rp(i?^,) 3 . The Heaviside function 9(y) 
ensures that the elastic interaction between grains only 
exists if they are in contact. For the leftmost (rightmost) 



granule in the chain, only the second (first) term on the 
right hand side of Eq. © is present. Initially the gran- 
ules are placed along a line so that they just touch their 
neighbors in their equilibrium positions (no precompres- 
sion), and all but the leftmost particle are at rest. The 
initial velocity of the leftmost particle (k = 1) is V\ (the 
impulse). We define the dimensionless quantity 



Q 



MiV 2 



a{R[) 



I \n+l/2 



(4) 



and the rescaled quantities Xk, t, mk, and R k via the 
relations 



^1 „,!/«■/ 



—a 1/n t 



y k = R[a 1/n x k , 

R' k = R[R k , M k = M x m k . 

Equation ([3]) can then be rewritten as 

m k x k = r fc „i(x fe _i - Xk) n ~ 1 0(x k -i - x k ) 
-r k {x k - x k+ i) n ~ 1 9{x k - Xfe+i), 



(5) 



(0) 



where a dot denotes a derivative with respect to i, and 



rk 



2R k R 



fe+i 



Rk + Rk+i 



1/2 



(7) 



The rescaled initial velocity is unity i.e., vi(t = 0) = I . 
The velocity of the k-th granule in unsealed variables is 
simply V\ times its velocity in the scaled variables. 



III. DECORATED CHAINS AND THE 
EFFECTIVE DYNAMICS 

There are several ways to construct a decorated chain. 
We shall consider decorated chains where small granules 
(all of the same size) alternate with larger granules. Our 
goal is to arrive at an analytically manageable effective 
description of a decorated chain. To be useful, such a de- 
scription must also be accurate in that it should capture 
(hopefully quantitatively) certain properties of the pulse 
dynamics. In particular, we attempt to capture the time 
dependence of the velocity profiles of the granules as the 
pulse moves along the chain, and the time that it takes a 
pulse to travel along the chain. We also try to character- 
ize the frequency of oscillation of the small granules that 
separate each pair of large granules. 

The frequency of oscillation of the small granules can 
be calculated by considering just a threesome of grains, 
two large ones with a small one in between. The anal- 
ysis of a three-grain chain is also a convenient way to 
detail the mathematical analysis that is then repeated 
for a longer system. The three-grain chain is analyzed 
in Sec. IIII Al On the other hand, to characterize the 
large granules that will constitute the elements of our 
effective chain, the shortest unit that must be analyzed 
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is one of five granules in the sequence large-small-largc- 
small-large. This is the minimal chain that starts with a 
large granule (as do our long chains), and has an interior 
granule with small granules on either side (whose effects 
will be captured in the renormalized mass and renormal- 
ized interactions). The five-grain chain is analyzed in 
Sec. IIII Bl Our effective chain of arbitrary length will 
then consist of all large granules, two end ones identi- 
fied with the two end granules in the five-grain sequence, 
and any number of interior granules identified with the 
interior large granule in the five-granule sequence. 

A. Three-grain chain 

We first consider a chain of three granules labeled k—l, 
k and k + 1, granule k being the small one. As noted 
above, this analysis will tell us something about the fre- 
quency of oscillation of the small granule, and it will also 
allow us to lay out the approximation scheme to be used 
in longer chains. The dynamics of the three granules is 
governed by the equations of motion 

Xk-i = -Tl(x k -i - x k ) n ' 1 B{x k -x - x k ) (8a) 

x k +i = U(x k - x k+ i) n ~ 1 6(x k - Xk+i) (8b) 
mx k = Tl(x k -i - Xk) n ~ X 6{xk-i - x k ) 

-K{x k -x k+1 ) n - x 0{x k -x k+1 ). (8c) 

To convey the method most clearly, we have here taken 
the masses m k ±i of the larger particles to be equal to 
each other and equal to unity in rescaled variables; m 
is the mass of the small granule. The effective elastic 
constant 1Z is given by Eq. |(7J) for a large granule of unit 
radius adjacent to a small one of radius r, that is, 

K =(— ) ' ' 9 > 

The generalization to tapered chains is implemented 
later. 

When granule fc — 1 is given an initial impulse, the 
energy is transported to granule k, which in turn transfers 
it to the last granule fc+1. Since granule k is small, it gets 
compressed between the two large granules and oscillates 
between them during the energy transfer process. The 
frequency of these oscillations can be expressed in terms 
of the maximum force that the granule k experiences from 
granule k + 1 or k — 1. We follow Melo et al. [35[ and 
write, 

x k {t) = Xk(t) + Asin(ut + <f>) (10) 

Here A, tu and (f> depend on the size ratio between the 
small and large granules and may also depend on the 
time. We assume that these time dependences are negli- 
gibly weak during the time of energy transfer from k—l 
to fc+1. The first term x k represents the average motion 
around which the particle executes oscillations of small 



amplitude A with frequency ui. Substituting Eq. (JTDJ) in 
Eq. (|8cp . and assuming that the amplitude A is much 
smaller than x k — x k +i and x k -i — x k , we expand in A 
and from the lowest two orders we obtain 

mx k w TZ(x k -i - x k ) n ~ 1 9{x k ^i - x k ) 

-Tl(x k -x k+ - l ) n ~ 1 6{x k -x k+ -s) (11a) 

muj 2 w TZ(n - 1) [(x k -i - x k ) n ~ 2 6(x k -i - x k ) 

+ {x k - x k+1 ) n - 2 6(x k - Xh +i)] ■ (Hb) 

In the following, for economy of notation we do not the 
write ^-functions explicitly, but their presence restricting 
the interactions to compressions should be kept in mind. 

The small granule oscillates with maximum frequency 
around the "equilibrium" position, defined as the posi- 
tion at which the forces F obtained using the average 
displacement of the smaller granule, cancel each other on 
the two sides, i.e., where x k = 0. Setting 

F = n{x k ^ - x,)™- 1 = TZ(x k - Xk+.Y 1 - 1 (12) 

[cf. Eq. (jllap j. and using Eq. (|llb[) . we can write for the 
maximum frequency / = cj/27t, 

1 l2(n - 1) i - ,»-2 

This is the theoretical prediction for the frequency of the 
small granule that we will later compare with numeri- 
cal simulation results once we are able to calculate the 
value of the average force F using our binary collision 
approximation. 

At the time when the small particle executes oscilla- 
tions with maximum frequency, we have x k = 0. From 
Eq. (|llap . this implies that 

x k = ^(x k +i + x k -i). (14) 
Also, note that from Eqs. Ij5a ]) -(|5c ]) 

mx k + Xk-i + x k +i = 0. (15) 

Equation (fl~5|) is the statement that the net force on 
the system is zero. Substituting Eq. (fTT)|) in (TTB]) . us- 
ing Eq. (|14p . and ignoring the small oscillatory part, we 
get (y + l) (xk-i +x k +i) = 0, which describes a system 
of two granules each with a renormalized mass 

M = l + |- (16) 

Replacing x k in Eqs. (|5a|) and (|8b[) by the form assumed 
in Eq. (TTUJ) and using Eq. (JT5J) we obtain 

/ \n— 1 

mA 2 . 
H tu sin + (p) 

i \n— 1 

-— — u 2 sii\{u)t + <j>). (17) 
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Since the oscillation amplitude A is small, wc neglect the 
oscillatory terms in Eq. (JTTJ) and simply write 



n 



fJ.Xk-1 



(J-Xk+1 



n-l 



2 



Xk+l 



,n-l 



(18) 



Equations (fT5)) define the effective interaction potential 
V e ff between granules k — 1 and fc + 1. Thus we can 
treat the system of two identical granules decorated with 
a smaller granule as an effective system of two granules 
(now relabeled as k and k + 1) with effective masses [i 
interacting with an effective potential given by 



eff 



n2 



K 



x k+1 ) n . 



(19) 



Thus the overall effect of an intermediate granule is to 
reduce the interaction strength between granules k — 1 
and k+l (in the original chain) and increase their masses. 
This shows that the pulse will be wider as compared to 
the monodisperse case. The result (|19[) remains valid as 
long as the oscillation amplitude of the small granule is 
small enough so that only the lowest order terms in the 
expansion, Eqs. pia[) and (|llb[) . are important. We shall 
see that this puts an upper limit on the size of the smaller 
granule. 



B. Long decorated chain 

Using the three-grain chain to study a tapered chain 
leads to a certain ambiguity because when there are two 
large granules of distinct masses around a small one, the 
renormalized masses will be different depending on the 
pairing. Thus, the renormalized mass of the fcth granule 
would have different values when considering the interac- 
tion between granules k — 1 and k or between k and k+l. 
This ambiguity can be suppressed if we consider a chain 
of five granules centered on granule k, labeled from k — 2 
to k + 2. Granules k — 1 and k + l are small ones of ra- 
dius r. The radius of large granule i = fc, fc± 2 is Ri> r. 
Moreover, this approach is useful even for a uniform dec- 
orated chain because it naturally shows the distinction 
between a granule inside the chain and a granule on ci- 
ther border. The dynamics of this chain of granules is 
governed by the set of equations 



m k - 2 Xk-2 = -r k -2(xk-2 - Xk-i 

\n-l. 



,n-l. 



Xk-2 ~ Xk-l) 

rriiXi = ri-i(xi-i - x l ) n ~ 1 9(x i ^i - x 4 ) 

- ri(xi — Xi + i) n ~ l 6(xi - £j + i), i = k, k±l, 
m k+2 x k+2 = r k +i(xk+i - x k+2 ) n ~ 1 0(xk+i ~ x k+2 )- 

(20) 

For the smaller granules, as was done in Eq. (p~0|) . we 
again assume a separation into an average displacement 
and an oscillatory contribution about the average. Fol- 
lowing the same steps that led from Eq. (JTUJ) to Eq. (fT4")) . 



we now obtain for the smaller granules 

Xi-l+U t X l+ l 



1 + «! 



where i = k — 1, fc + 1, and 



(21) 



(22) 



Using Eq. ([21]) in Eq. (|20|) . the pulse propagation in 
a decorated chain of five granules can be described by 
pulse propagation along an effective chain of three large 
granules, the "left" (I), the "middle" (m), and the "right" 
(r), of renormalized masses 



/'/ 



m k 



m k -2 + 



rk-2 
Ho-l 



rk+l 
r k 



Tk-l 

rk-2 



m k+2 + 



1 



'•fc 
rk+l 



(23a) 
(23b) 

(23c) 



Henceforth we relabel the three grains in the effective 
chain as (I, m, r) — > (k — 1, fc, k + 1). The effective inter- 
action between the fcth and (k + l)st granules (relabeled) 
is given by 



Ve// = — (a*-a*+i) n , 



where 



C*(n) = 



1 



(24) 



(25) 



and where r k is now given by 



r k 



2R k r 
Rk+r 



(26) 



A long decorated chain can now be represented by an 
effective chain of large granules where the mass of the 
fcth granule (except for those at the edges of the chain) 
is fx m . Note that the effective mass /x m is modified from 
m k by the presence of the two small granules on cither 
side of it in the decorated chain. The masses of the gran- 
ules at the edges of the long effective chain are /i; and 
H r and are modified from m k ^2 by only a single small 
mass. The interaction between the fcth and (fc+ l)st (re- 
labeled) granules in the effective long chain is given by 
Eq. (j2~4")) with (|25[) . Hence, as in the three-grain chain, 
the interactions are weaker and the masses larger in the 
effective chains than in the original decorated chain. For 
the special case, R k = R k ±2 in the five-granule decorated 
chain, Eqs. (HU) and (gj) reduce to Eqs. (JUJ) and 1(15]). 
respectively. 
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IV. BINARY COLLISION APPROXIMATION 
FOR EFFECTIVE CHAIN 

We now use the binary collision approximation to ob- 
tain analytic expressions for the pulse propagating along 
a decorated chain and make the connection with the ef- 
fective description of Sec. IIIII The following analysis is 
based on Ref . [32j . 

In the binary collision approximation we assume that 
the transfer of energy along the chain occurs via a succes- 
sion of two-particle collisions. First, particle fc = 1 with 
velocity V\ — 1 collides with initially stationary particle 
k = 2, which then acquires a velocity vi and collides with 
stationary particle k = 3, and so on. Using energy and 
momentum conservation and remembering that in the ef- 
fective chain the granules have effective mass fj,k, it can 
easily be shown that the velocity Vk of the fcth granule is 
given by 



fc-i 



n.r 



fiy 



(27) 



The coupled dynamics of any two interacting granules 
in our models can be reduced to a single particle dynam- 
ics. This is helpful toward obtaining analytical expres- 
sions for some of the pulse properties, as we shall see 
below. For this purpose, we pick granules k and k + 1 
and introduce the difference variable 



Zk = Xk - Xk+i- 



(28) 



The equation of motion for the difference variable is ob- 
tained by subtracting the equations of motion of the two 
granules during a collision, 



Xk 



Xk+l 



ftk+l 



(X k - Xk+lY 



(x k - Xk+lY 



(29) 



where Cfc( n ) is defined in Eq. 
(|28[) immediately leads to 



where 



Zk 



Mi 



Ck(n) 



Equation (|2SJ) with 



(30) 



fi-k + 



(31) 



is the reduced mass of the coupled system. 

Equation (|3"0|) is the equation of motion of a particle of 
effective mass Aik in the potential C,k{n)z^/n defined for 
Zk > 0. The initial conditions are ifc(0) = Vk since the 
velocity of granule k + 1 is zero before the collision, and 
Zfc(0) = since there is no precompression. The energy 
conservation condition 



Ait) 



nMk 



Zk(t) 



1 



(32) 



leads to 



Zk{t) 



£(o) 



Kk{n) 



1/2 



(33) 



We say that the pulse arrives at granule k when the 
velocity of granule k surpasses that of granule k — 1, and 
that it moves on to granule k + 1 when the velocity of the 
(k + l)st granule surpasses that of the fcth granule. The 
residence time Tk on granule fc is the time that granule 
fc takes to transfer the pulse from fc — 1 to fc + 1, and is 
given by 



T k = 



dz k _ 
z k 



dz k 



1/2 ; 



(34) 

where z™ ax is the compression when the velocities of par- 
ticles fc and fc + 1 are equal (which is also the maximum 
compression) , 



nMk . 2 



(0) 



l/n 



(35) 



The integral can be performed exactly to yield 



nMi 



l/ra 



OA ( \ I ^fc(O) 1 — 7-7-r. 36) 



Finally, the total time taken by the pulse to pass the fcth 
granule is obtained by summing Tk, 



(37) 



In the next section we compare these various quantities 
obtained from the binary collision approximation imple- 
mented on the effective undecorated chains with deco- 
rated chain numerical integration results. 

At this point we can also calculate the average force 
that appears in Eq. (TT3"]) . When the two granules fc and 
fc + 1 collide, the maximum force between them corre- 
sponds to the maximum compression z™ ax . Thus the 
maximum force is 



f = aw^rT" 1 

nM fc if(0) 



l/r, 



(38) 



The explicit formulas for the three-grain configuration 
considered in Sec. IIII Al correspond to using Eq. (f3i"j) with 
fi k = Hk+i = ft and hence M k = ft/2, Cfc(n) = TZ/2 n ~ 1 , 
and i^(0) = 1. We obtain 



F= VfL 

4 Vn2 



7^ 



l/n 



(39) 
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where 1Z and /i are given in Eqs. © and (flB)) . respec- 
tively. Substituting this in Eq. (|13p . we obtain 



/ 



J_ 2{n-l)n fll 
2ir V m \ /i 



(40) 



In the next section we compare this frequency with that 
obtained from the numerical solution of the exact dy- 
namical equations (|5a|) - (|8c|) . 



V. COMPARISON WITH NUMERICAL 
RESULTS 

While our binary collision approximation is valid for 
n > 2 [l2j], in our numerical comparisons we focus on 
spherical granules. We start with the chain of three par- 
ticles of Sec. IIII Al Our first figure, Fig. [1] shows that 
the important step of neglecting oscillatory terms as de- 
tailed in the derivation of Eqs. (|lla[) . (|llb|) and (fT5|) is 
justified. The figure shows the displacement of the three 
granules as a function of time. The dashed curve shows 
the average displacement, Eq. (fl"4|) . of the small gran- 
ule. The "small amplitude" in question is the difference 
between the oscillatory curve and the dashed curve com- 
pared to the differences between the dashed curve and 
the two large-granule displacement curves. 




FIG. 1: (Color online) Solid lines: displacements of the 
three granules in a three-granule chain as a function of time. 
Dashed line: average displacement of the small middle gran- 
ule. 

Figure [5] shows the velocity profile of the granules ob- 
tained from the solution of the exact dynamical equa- 
tions for the three granules. The small granule oscillates 
a number of times before the chain disintegrates. In the 
same figure we also show a comparison with the results 
obtained from the effective dynamical equations (|18| . As 
we see from the figure, the behavior of the large particles 
is captured very accurately by the effective interaction 
given by Eq. (flT)|) . We observe this quality of agreement 
as long as the radius of the small granule is less than 
about 40% of the larger ones. We return to this point 





1 ? 




1.0 




U.o 


>city 


U.o 


<j 


0.4 


> 


0.2 




0.0- 




-0.2 









1 2 3 4 5 

t 

FIG. 2: (Color online) Typical velocity profiles of three gran- 
ules (continuous curves) obtained from the numerical solution 
of the exact dynamical equations. The small particle of ra- 
dius r = 0.3 oscillates with almost constant frequency while 
the amplitude of oscillations shows a weak time dependence. 
The dashed curves represent results from the effective dynam- 
ics. 



later in the paper. The effective description improves as 
the size of the small granule decreases. 

Although it is not necessary for the characterization of 
the propagation of the pulse, it is interesting to compare 
the theoretical [Eq. (|40|) ] and numerical results for the 
oscillation frequency of the small particle as a function of 
its radius. The theoretical result is shown as a solid curve 
in Fig. [3] The numerical results are obtained by fitting 
the derivative of Eq. flO]), v(t) = 0.5 + u;Acos(ujt + </>)) 
(the particle velocity oscillates around 0.5). The resulting 
frequencies are shown as dots in the figure. While the 
theoretical curve appears closer to the numerical results 
for the larger granules, the percent difference between 
the two is 14% for the smallest granules and 20% for the 
largest granules shown in the figure. 
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FIG. 3: (Color online) The dots represent the numerical re- 
sults for the average angular frequency of oscillations of the 
small particle during the time of energy transfer from granule 
k — 1 to k -+- 1. The solid curve is the theoretical estimate 
based on the binary collision approximation using the effec- 
tive interaction. 
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A. A simple decorated chain 

A simple decorated chain is one in which the large 
granules all have the same radius. For our numerical sim- 
ulations of a simple decorated chain we consider TV = 50 
large granules, with each pair of large granules separated 
by a small granule. Initially all the granules are placed 
such that they barely touch their nearest neighbors (no 
precompression) . In all decorated chains, we consider the 
first granule to be a large one. A schematic of a typical 
simple decorated chain is shown in Fig. |U 




FIG. 4: Schematic of a simple decorated chain of nine gran- 
ules. 

Figure [5] shows the velocity pulse propagating through 
the chain. The smooth higher velocity peaks are associ- 
ated with the large granules, while the lower peaks with 
oscillatory contributions are those of the small granules. 
The oscillations in the velocity profiles of the small parti- 
cles are most pronounced around their maximum veloc- 
ity. This is due to the fact that the forces on the two sides 
of the small granule tend to cancel each other around 
that time. Almost all the noisy behavior at the bottom 




t 



FIG. 5: (Color online) Velocity profile of granules in a simple 
decorated chain with small granules of radii r = 0.3. 

of Fig [5] is due to the oscillations of the smaller granules. 
None of these oscillatory and noisy contributions affect 
the smoothness of the pulse front located on the larger 
granules. We also note that the initial energy given to 
the left edge of the chain is transported mainly through 
the larger granules. The energy of the small granules is 
negligible compared to that of the larger ones. 

In Fig. [6] we compare the exact dynamics of the larger 
granules with the results obtained from the solution of 
the dynamical equations of the associated effective un- 
decorated chain. We show only a small stretch in the 



figure, but the agreement is typical of the entire chain. 
A schematic of this chain is shown in Fig. [JJ The dynam- 
ical equations for the two edges of the effective chain are 
given in Eq. (JTSJ) with the associated effective mass (|16p , 
or cquivalcntly, the effective mass (I23bl) or (I23cl) for equal 
large granules. The interior masses satisfy the dynamical 
equation (with ^-functions omitted) 

Xk = 2^ K**-i - a*)"" 1 - (** ~ ^-i)' 1 " 1 ] (41) 

with the effective mass fl23ap , which for the simple deco- 
rated chain reduces to 

H k = n = 1 + m. (42) 

For clarity, we only show a section of the chain between 
k — 8 and k — 24 (here k counts only the larger granules). 
The results from the effective dynamics are in excellent 
agreement with the exact solution. 




t 



FIG. 6: (Color online) Exact results (solid curves) for the 
velocity profiles of the larger granules in a simple decorated 
chain for r = 0.3. The dotted curves show results from the 
solution of the effective dynamical equations. 




FIG. 7: Schematic of the effective undecorated chain asso- 
ciated with Fig. [4] Each pair of granules interacts with the 
same pair potential V e f / but the end masses are smaller than 
the interior masses. 

For a simple decorated chain, Tk is independent of k 
(except for the edge granules, which have a different ef- 
fective mass) and is given by 

/ n2"- 3 M \ 1/w r(l + l/») 

T k = ^[-n^r) r(i/n + i/2)' (43) 
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with the effective mass of Eq. (|42j) . The time t defined 
in Eq. (f57j) therefore varies linearly with k. The slope of 
the line is determine by T^. In Fig. [8] we show numerical 
results for the time taken by the pulse to reach the fcth 
granule. Note that k here labels only the larger particles 
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FIG. 8: Time t taken by the pulse to reach the fcth granule 
in a simple decorated chain for r = 0.1 (open circles), r = 0.2 
(filled circles), r — 0.3 (squares), and r = 0.4 (plus signs). The 
theoretical results are shown with dashed (r = 0.1), dashed- 
dotted (r = 0.2), dotted (r = 0.3), and continuous (r = 0.4) 
lines. 

in the actual chain, and for clarity we show results for 
the restricted range 20 < k < 50. Various symbols rep- 
resent different values of r. In the same figure we also 
show the analytical result Eq. ([57)) of the binary collision 
approximation. No fitting parameter are involved, and 
the agreement of the two results is again gratifying. This 
shows that the effective binary theory can give quantita- 
tively reliable results for simple decorated chains. 

These quantitatively reliable results are obtained pro- 
vided the radius of the small granules is no larger than 
about r ~ 0.4. In Fig. [5] we show the percent difference 
in the pulse amplitude calculated from the decorated dy- 
namics and from the effective undecorated one as a func- 
tion of the radius of the small granules. The percent 
difference is calculated by taking the percent difference 
in the pulse amplitude for each granule in the chain, sum- 
ming over all granules, and dividing by the total number 
of granules. The average percent difference remains un- 
der 1% up to about r = 0.4 and beyond that rises sharply. 



B. A tapered decorated chain 

We next consider a linearly tapered chain decorated 
with small granules of the same size. For numerical pur- 
poses, we again consider a tapered chain of N = 50 
larger granules. The radius of the kth granule is Rk = 
1 — S(k ~ 1) and S is the tapering parameter 32]. This 
chain is then decorated with the smaller granules of ra- 
dius r. A schematic of a typical tapered decorated chain 
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FIG. 9: Percent difference between the exact dynamics and 
the effective dynamics as a function of the radius of the small 
granules. The percent difference is an average calculated by 
taking the percent difference in the pulse amplitude for each 
granule in the chain, summing over all granules, and dividing 
by the total number of granules. 



is shown in Fig. 1101 




FIG. 10: Schematic of a tapered decorated chain of nine gran- 
ules. 



Our effective description is valid as long as the am- 
plitudes of oscillations of the smaller particles are small 
compared to their average displacement. This again oc- 
curs as long as the size of the small granules is roughly 
less than 40% of the size of the larger granules. The value 
of S is therefore restricted by the length of the tapered 
chain. This is the reason for choosin g lin ear rather than 
geometric tapering for this analysis [32j, which is even 
more restrictive. 

In Fig. [IT] we show the change in the pulse amplitude 
with the granule number. As already noted in the case of 
undecorated tapered chains [12] , the absolute value of the 
pulse amplitude is not captured correctly by the binary 
approximation, but the rate of change of the amplitude 
is close to that obtained from the exact results. 

However, other characteristics of the pulse are ex- 
tremely well captured by the binary approximation. We 
next compare the binary prediction for the time spent 
by each granule in the pulse. The analytical result for 
the residence time of the pulse on granule k is given by 
Tfe in Eq. (f3~4")l . In Fig. [J2] we compare the results ob- 
tained from the exact dynamics to those of the binary 
collision approximation for r = 0.1, 0.2, and 0.3. Ex- 
cellent agreement is again obtained between the two re- 
sults. The small initial disagreement is due to an edge 
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FIG. 11: Change in the pulse amplitude as a function of k. 
The filled and the empty circles represent the exact numerical 
result and the binary collision approximation result, respec- 
tively. The three different data sets are for S = 0.002, 0.004, 
0.006 from bottom to top. In all cases r — 0.3. 



effect. The residence time of the pulse decreases rapidly 
as the pulse propagates through the chain, as it does in a 
non-decorated forward tapered chain [HI ■ Note that the 
residence time of the pulse on granule k decreases with 
increasing r because the increase in the effective mass is 
accompanied by an even greater increase in the effective 
interaction bewteen granules. In the inset of Fig. we 
show the results for r = 0.4. Here the exact and approx- 
imate results disagree for large k values because the size 
difference between the larger and the smaller granules de- 
creases to where the description in terms of an effective 
chain breaks down. 



VI. CONCLUSIONS 



We have studied pulse propagation in decorated chains, 
that is, chains in which small and large granules alter- 
nate. We have presented a method to map the actual 
chain into an effective undecorated chain. This effective 
description allows us to obtain analytical results using the 
binary collision approximation for quantities that char- 
acterize the pulse properties in decorated chains. The 
effective description is extremely accurate as long as the 
size ratio of the small granules to the large granules is less 
than ~ 0.40. We have explicitly considered two different 
chains, simple decorated chains and forward tapered dec- 
orated chains. For simplicity, in order to obtain an ef- 
fective description, we have assumed a fixed size for the 
small granules. However the method can be easily ap- 
plied with little modifications to other more complicated 
decorated chains where small granules are also allowed to 
change size. The method is also implcmcntable to deco- 
rated chains with more than one large granule between 
the small ones. A more extensive generalization is nec- 
essary to deal with decorated chains with more than one 
small granule between large ones and thence to other size 
profiles. This work is in progress. 
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FIG. 12: Residence time Tk of the pulse on the fcth granule 
for r — 0.1, 0.2, and 0.3 from top to bottom, and S = 0.01. 
Numerical and analytical results are represented by the filled 
and open circles, respectively. The inset shows results for 
r = 0.4. 
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